F/6  9/2 


AD-AO06  329 
UNCLASSIFIED 


ROYAL  AIRCRAFT  ESTABLISHMENT  FARNBOROU6H  (ENGLAND) 

A  STUDY  OF  LANDSAT  MULT I -SPECTRAL  SCANNER  DATA  TAPES. (U) 


FEB  80  A  H  BENNY 
RAE-TR-8001B 


NL 


ORIC-BR-73767 


adA0863 


ROYAL  AIRCRAFT  ESTABLISHMENT 


Technical  Report  80018 
February  1980 


A  STUDY  OF  LANDSAT 
MULTI-SPECTRAL  SCANNER 


DATA  TAPES 


A.H.  Benny 


DTIC 

ELECTE 
JUN2  5 


Procurement  Executive,  Ministry  of  Defence 
Fernborough,  Hants 


UNLIMITED 

80  6  19  1 


'  SUMMARY 

A  detailed  study  is  described  of  a  number  of  computer  compatible  magnetic 
tapes  containing  Landsat  multi-spectral  scanner  (MSS)  data.  This  study  reveals 
the  differences  between  raw  and  corrected  data  tapes  and  indicates  how  the  con¬ 
version  is  performed.  It  also  shows  the  differences  between  the  existing  and  a 
new  form  of  corrected  data  tape. 

From  the  knowledge  gained,  a  computer  program  has  been  written,  to  enable 
raw  data,  desirable  for  certain  computer  analyses,  to  be  obtained  from  the 
existing  stocks  of  corrected  data  tapes. 
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INTRODUCTION 


Amongst  the  many  satellites  at  present  orbiting  the  Earth  is  a  series  known 
as  Landsat.  These  vehicles  contain  a  number  of  sensors  which  make  observations 
of  the  surface  of  the  Earth,  the  information  obtained  being  broadcast  by  the 
satellites  and  collected  at  ground  receiving  stations.  The  type  of  sensor  of 
interest  in  the  present  study  is  the  multi-spectral  scanner  (MSS),  which  is 
described  in  Ref  1 . 

Normally,  the  received  data  is  recorded  on  high  density  magnetic  tapes 
which  are  not  computer-compatible,  that  is  they  are  not  capable  of  being  read 
directly  into  a  conventional  computer  installation.  It  is  therefore  common  for 
ground  receiving  stations  to  have  the  special  equipment  necessary  to  convert  high 
density  tapes  to  the  more  widely  usable  computer  compatible  tapes  (CCTs). 

The  conversion  from  high  density  tape  to  CCT  may  be  done  with  little  or  no 
alteration  to  the  data  content,  and  such  a  CCT  may  be  thought  of  sub  containing 
•raw'  or  unprocessed  data.  Since  most  users  require  corrected  data  tapes,  it  is 
common  practice  for  the  high  density  tapes  to  be  converted  directly  into 
•corrected'  CCTs,  without  intermediate  storage  in  'raw'  form. 

Receiving  stations  may  therefore  offer  to  supply  Landsat  MSS  data  users 
with  either  raw  or  corrected  data  tapes.  Until  recently,  the  demand  has  been 
solely  for  corrected  tapes.  For  exsimple,  Space  Department ,  RAE  has  acquired  a 
considerable  number  (about  200)  of  corrected  CCTs  but  had  no  raw  CCTs  at  the  time 
that  this  work  was  started. 

During  the  course  of  computer  analysis  of  Landsat  MSS  data  it  has  become 
apparent  that  corrected  data  suffers  from  certain  adverse  features,  resulting 
from  the  correction  process.  It  was  therefore  considered  desirable  to  obtain  raw 
data,  and  in  view  of  the  considerable  amount  of  corrected  data  held  at  RAE  it 
seemed  worthwhile  to  attempt  to  devise  a  method  (ie  a  computer  program)  to 
restore  raw  data  from  corrected  tapes.  This  Report  describes  a  study  which  has 
been  made  to  determine  the  data-structure  of  raw  and  corrected  tapes  and  the 
relationship  between  them,  with  the  aim  of  converting  corrected  data  back  to  the 
raw  form. 

1.1  Cautionary  note 

Many  scientific  investigations  are  concerned  with  natural  objects  or 
phenomena,  whereas  this  study  is  concerned  with  man-made  objects  (CCTs).  There 
is  a  possibility  that  a  purely  descriptive  observation  nay  be  interpreted  as 
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implying  adverse  criticism.  The  author  wishes  to  make  it  clear  that  this  Report 
is  not  intended  in  any  way  as  a  criticism  of  the  tapes  under  study. 

2  DESCRIPTION  OF  COMPUTER  COMPATIBLE  TAPES 

The  'corrected*  computer  compatible  tapes  are  described  in  some  detail  in 
Ref  2,  dated  May  1976  and  redefined  in  Ref  3  dated  October  1979»  which  became 
available  during  the  course  of  the  work  described.  These  publications  give  con¬ 
siderable  detail  about  the  tapes,  though  they  leave  a  number  of  small  but  impor¬ 
tant  points  unresolved.  No  description  of  'raw'  data  is  given,  and  insufficient 
details  of  the  conversion  from  raw  to  corrected  data  are  provided  to  allow  the 
process  to  be  reversed. 

A  comparison  of  Refs  2  and  3  shows  that  significant  differences  exist 
between  the  relevant  formats,  so  this  study  will  refer  to  'old'  and  'new'  correc¬ 
ted  tapes,  the  former  being  produced  before  the  summer  of  1979  and  the  latter 
after  that  time.  On  the  basis  of  information  received  by  private  communication, 
the  changeover  date  has  been  selected  as  1  August  1979- 

When  the  present  study  began,  it  was  not  known  that  a  format  change  was 
about  to  occur.  It  is  perhaps  fortunate  that  the  change  occurred  at  this  time, 
as  this  tudy  was  then  extended  to  cover  the  new  format,  and  has  been  able  to  be 
more  complete  in  itself  and  in  the  output  products  (ie  computer  programs)  which 
have  resulted  from  it. 

2.1  General  structure  of  the  CCTs 

The  general  structure  of  all  CCTs,  both  raw,  and  old  and  new  corrected,  is 
similar,  and  is  illustrated  in  Table  1.  The  magnetic  tapes  contain  a  sequence  of 
'records',  as  shown,  starting  with  two  header  records,  separated  by  a  tape 
marker.  The  next  six  records  refer  to  transformations  and  calibrations  relevant 
to  the  scene,  and  are  followed  by  a  tape  marker.  There  are  9144  data  records, 
referring  to  the  four  spectral  bands  of  data  for  the  2286  lines  or  scans  of  the 
scene. 

Refs  2  and  3  describe  in  some  detail  the  nature  and  contents  of  these 
records.  This  Report  is  concerned  only  with  details  of  the  records  not  clarified 
by  these  references,  and  with  differences  between  the  various  raw  and  corrected 
tapes.  If  it  is  intended  to  study  section  3  in  detail,  it  is  advisable  that 
either  Ref  2  or  3  be  available  for  concurrent  use. 


2.2  The  data  tapes  examined 


Space  Department,  RAE  had  an  extensive  collection  of  corrected  CCTs  all  of 
the  'old'  type,  when  this  study  started.  Two  more  tapes  were  obtained  especially 
for  this  study,  these  being  a  raw  and  a  corrected  tape  (nominally  of  the  same 
scene).  As  the  study  proceeded,  it  became  apparent  that  the  corrected  tape  was 
of  a  different  and  new  type.  Since  then,  several  other  tapes  have  been  examined, 
the  details  being  shown  in  Table  2.  All  of  these  tapes  were  processed  at  Fucino 
in  Italy,  over  a  period  of  some  21  months.  The  tapes  referred  to  as  D  and  E 
contain  data  from  Landsat  1,  the  others  from  Landsat  2.  The  scene  numbers  refer 
to  the  satellite  path  and  row  .  Tapes  A  and  B  are  the  same  scene,  a  portion  of 
Sweden,  and  both  derive  from  the  same  satellite  pass.  The  other  scenes  are  of 
the  south-west  of  England,  including  the  Severn  estuary. 

3  DETATT.TO  INVESTIGATION  OF  THE  TAPES 

The  computer  compatible  magnetic  tapes  were  examined  in  detail,  on  a 
Prime  400  computer  system  including  a  1600  bpi  9-track  magnetic  tape  unit. 

A  computer  program  was  available  which  allowed  the  tapes  to  be  studied  in 
close  detail.  However,  due  to  the  very  large  quantity  of  data,  it  was  not  pos¬ 
sible  merely  to  print  it  all  out  and  study  the  listing,  as  this  would  have  run 
to  many  thousands  of  closely  packed  pages  for  each  tape.  Instead,  severed  com¬ 
puter  programs  were  written,  to  explore  various  facets  of  the  information 
required.  Some  of  these  programs  are  mentioned  in  the  course  of  this  and  the 
next  section. 

This  section  is  not  to  be  taken  as  a  complete  description  of  the  tape 
formats,  for  which  attention  should  be  turned  to  Ref  2  or  3*  Mention  will  only  be 
made  here  of  detail  not  noted  in  those  references. 

3.1  JSC  header  record 

The  JSC  or  Johnson  Space  Centre  Header  Record  contains  3060  bytes,  most  of 
which  are  unused.  Only  two  items  have  been  found  in  the  JSC  header  on  all  tapes, 
but  some  tapes  have  a  few  additional  items. 

The  items  specified  in  Ref  3  a re  listed  in  Table  3»  All  of  the  tapes 
studied  contain  items  1  and  2.  It  is  of  interest  that  the  values  provided  on 
tapes  A  and  B,  which  are  nominally  the  same  scene,  differ  slightly.  This  may  be 
because  the  scene  is  in  fact  slightly  displaced,  between  tapes  A  and  B,  or  it  may 
reflect  a  more  accurate  value  in  the  case  of  the  corrected  tape  B. 


6 


The  old  corrected  tapes  examined  sill  contain  item  3  and  the  new  corrected 
tapes  also  contain  item  4. 

3.2  Landsat  header 

Whilst  the  JSC  header  always  uses  EBCDIC  characters,  the  Landsat  header 
and  subsequent  records  may  use  either  EBCDIC  or  ASCII.  The  character  code  is 
indicated  by  one  of  the  process  flag  bits  in  the  Landsat  header.  All  of  the 
tapes  examined  have  been  found  to  use  EBCDIC  characters. 

The  contents  of  the  Landsat  header  are  described  in  Refs  2  and  3 •  The 
differences  noted  between  the  headers  on  the  raw  tape  (A)  and  corrected  tape  (B) 
for  a  nominally  identical  scene  are: 

i 

Value 


Item 

Contents 

Raw 

Corrected 

1 

Day  number  since  launch 

899 

1899 

2 

Frame  identification 

09012 

2208018049 

3 

Centre  latitude  N 

5990 

5984 

4 

Centre  longitude  E 

-1799 

1768 

5 

Cycle  number 

50 

49 

6 

Date  master  tape  generated 

5  April  79 

16  Sept  79 

7 

Date  copy  tape  produced 

5  April  79 

16  Sept  79 

8 

Process  flags 

0000000 

1111010 

Comments  may  be  made,  on  the  differences  between  the  two. 

Item  1  Raw  data  omits  leading  *  1  * :  this  is  the  true  value,  whereas  the 
corrected  value  is  in  error. 

Item  2  The  frame  identification  on  corrected  tapes  includes  the  mission 
(Landsat  2)  track  (208)  frame  (01 8)  and  cycle  (049). 

On  the  raw  tape  it  may  refer  to  hours  (9)  minutes  (01)  and  tens  of 
seconds  (2)  at  which  the  centre  point  was  imaged,  ie  9:01:20  on  the  day 
specified. 

Item  3  The  small  differences  here  probably  represent  the  slightly  different 
and  4  area  imaged  (see  section  4.1).  The  differences  are  of  the  required 

order  of  magnitude,  for  a  satellite  near  to  the  northerly  apex  of  its 
orbit,  assuming  that  the  values  represent  degrees  (two  figures)  and 
hundredths  of  degrees  (two  figures).  The  specifications  state  that  the 
second  pair  of  figures  represent  minutes  of  arc,  but  this  clearly  cannot 
be  the  case,  as  the  values  exceed  the  maximum  possible  value  of 
39  minutes. 


The  minus  sign  in  the  raw  data  for  item  4  is  not  understood:  it  is 
clearly  in  error. 

Item  5  It  is  not  known  why  these  values  differ.  The  actual  cycle  number  is  49, 
which  means  that  the  raw  data  tape  is  incorrect. 

Items  6  These  are  the  dates  at  which  the  CCTs  were  prepared,  and  of  course  need 
and  7  not  be  the  same. 

Item  8  Raw  data  is  ' unprocessed '  and  would  be  expected  to  have  process  flags 
all  zero.  The  code  for  the  processed  tape  implies:  (left  to  right). 

Bit  1 :  No  meaning  to  user 
Bit  2:  256  radiometric  levels 

Bit  3:  Geometric  correction  applied 
Bit  4:  No  meaning  to  user 

Bit  3'  Statistical  radiometric  correction  applied  (as  distinct  from 
step  wedge  calibration). 

Bit  6:  No  meaning  to  user 

Bit  7:  EBCDIC  character  code  (as  distinct  from  ASCII) 

It  is  noteworthy  that  the  other  corrected  tapes  have  the  same  process 
flags  as  tape  B. 

3.3  Transformation  record 

The  transformation  record  contains  up  to  36  values,  but  in  the  tapes 
studied  so  far  only  seven  items  have  non-zero  values.  The  values  on  the  raw  and 
corrected  tapes  are: 

Item  Meaning  Raw  Corrected 

1  UTM  zone  number  33  33 

2  UTM  Northing  (metres)  6643843.5  6636431.27 

3  UTM  Easting  (metres)  667227.25  650377.3597 

4  Frame  orientation  0.4261635542  0.275119713 

5  Pseudo  altitude  of  satellite  (km)  249.7637482  248.7237 

7  X  scale  factor  3.703999996  3.703999996 

8  Y  scale  factor  925.124939  3.703999996 

Items  2,  3  and  4  would  be  expected  to  differ  somewhat,  because  of  the 
difference  between  the  scene  framing,  mentioned  earlier.  Item  5  may,  or  may  not, 
depend  on  the  same  reason. 


Items  7  and  8  refer  to  the  70  mm  film  images  which  may  be  produced  as  an 
alternative  to  CCTs,  so  probably  have  no  relevance  to  the  computer  processing 
done  at  RAE. 

3.4  Radiometric  look-up  table  records 

Biere  are  five  look-up  table  records,  corresponding  to  the  spectral  bands 
4  to  8  inclusive.  This  study  is  concerned  only  with  bands  4  to  7,  so  band  8  is 
not  considered  here. 

The  look-up  tables  show  the  conversion  which  has  been  used  to  transform 
6-bit  values  (0  to  63)  from  the  satellite  sensors  to  8-bit  values  (0  to  255)  on 
the  corrected  tapes.  The  conversion  takes  account  of  the  known  behaviour  of  the 
sensors  (ie  bands  4  to  6  are  quasi-logarithmic,  whereas  band  7  is  linear)  and 
also  the  differences  in  sensitivity  between  each  of  the  six  detectors  associated 
with  each  spectral  band. 

For  any  scene,  a  set  of  look-up  tables  may  be  prepared.  The  6-bit  data  on 
the  high-density  tapes  is  either  used  unaltered  (for  raw  tapes),  or  transformed 
to  8-bit  data  by  means  of  the  look-up  tables  (for  corrected  tapes). 

The  contents  of  the  look-up  tables  have  been  examined  on  the  various  types 
of  tape,  and  the  differences  are  described  below. 

3.4.1  Raw  tapes 

For  raw  data  tape  A  all  of  the  look-up  tables  contain  only  zero  values. 

(This  is  to  be  expected  since  raw  data  does  not  need  such  tables.) 

3.4.2  Corrected  tapes 

For  corrected  tapes,  complete  look-up  tables  are  provided,  as  expected. 

Ref  2  explains  that  the  tables  are  derived  either  from  calibration  data  from  the 
satellite  (calibration  wedge),  or  from  a  statistical  analysis  of  the  scene  itself. 
Bit  5  of  the  process  flag  word  in  the  Landsat  header  indicates  that  the  latter 
has  been  done  for  all  corrected  tapes  studied. 

It  is  of  interest  to  note  that  the  look-up  tables  contain  a  few  instances 
of  two  adjacent  values  being  the  same.  Thus  it  is  not  always  possible  to  perform 
the  inverse  mapping  from  corrected  to  raw  values  without  an  occasional  ambiguity 
occurring.  This  does  not  occur  on  all  images,  and  it  is  possible  to  restore  data 
to  its  raw  form  in  many  cases,  at  least  as  far  as  the  radiometric  values  are 
concerned. 


3.5  Multi-spectral  scanner  data  records 

A  scene  is  defined  as  2286  lines.  Each  line  consists  of  four  spectral 
bands,  each  of  which  is  contained  in  one  magnetic  tape  record.  There  should 
therefore,  be  2286  x  4  =  91^  data  records  on  the  tape,  followed  by  a  file 
marker.  At  1600  bpi  this  almost  fills  a  2400  foot  tape. 

The  data  structure  may,  therefore,  be  considered  as  consisting  of  sets  of 
four  records  (bands  4,  5,  6  and  7)  for  each  line  of  the  scene.  Of  each  set  of 
four,  the  first  record  differs  from  the  other  three,  in  the  manner  noted  below. 

Each  record  consists  of  3780  bytes.  Bytes  1  and  2  form  a  16-bit  integer, 
having  the  values  1,  2,  3i  ^  for  each  of  the  records  of  a  set  of  four.  Thus, 
byte  1  of  the  record  is  always  zero,  and  byte  2  of  each  band  4  record,  for 
example,  would  have  the  value  1. 

For  band  4  data,  bytes  3-l80  contain  'ancillary  data'  (see  below)  and  bytes 
181-3780  contain  video  data.  For  the  other  three  spectral  bands,  the  video  data 
is  contained  in  bytes  3-3602  and  the  ancillary  data  is  in  bytes  3603-3780,  ie 
after  the  video  data.  The  ancillary  data  for  bands  5*  6  and  7,  where  it  exists 
is  usually  the  same  as  that  for  band  4. 

3.5.1  Ancillary  data 

The  raw  data  tape  A  contains  ancillary  data  only  in  the  band  4  record;  in 
the  other  three  bands  the  words  are  all  zeros. 

The  interpretation  of  the  ancillary  data  is  described  in  Refs  2  and  3*  but 
the  meaning  of  several  items  is  not  entirely  clear.  For  instance: 

'End  of  video  flag  count'  -  on  the  raw  data  tape  this  has  a  value  which 
may  be  interpreted  as  -19700. 

'Minor  frame  sync  loss  count'  -  zero  in  all  cases  examined. 

'Uncorrected  line  length'  -  this  seems  to  be  approximately  but  not  exactly 
equal  to  the  number  of  pixels  in  the  line  of  raw  data  which  corresponds  to  the 
corrected  line.  On  the  raw  data  tape  the  value  in  the  ancillary  block  is  zero. 

The  meaning  of  'start  pixel  number'  and  'end  pixel  number'  seems  to  vary 
from  one  type  of  tape  to  another.  This  is  described  in  more  detail  in 
section  3*5.2. 

Summarizing  the  ancillary  data  therefore: 

(a)  Raw  tapes  only  contain  ancillary  data  for  band  4;  for  other  bands  it 
is  all  zero.  Even  in  the  band  4  data,  only  a  few  of  the  values  are  present 
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namely:  time  and  day/time,  EOV  flag,  scan  number,  start  and  end  pixel  numbers, 
and  sensor  number. 

(b)  Old  corrected  tapes  contain  ancillary  blocks  for  all  bands,  but  those 
for  bands  5,  6  and  7  are  less  complete,  eg  they  lack  scan  line  numbers,  start  and 
end  pixel  numbers,  time  and  day/time.  For  band  4  the  EOV  flag  is  set  to  -19700 
(as  with  the  raw  data),  and  the  XSCAN  is  set  to  zero.  The  'one  fills'  are  also 
set  to  zero  but  most  of  the  other  ancillary  data  is  provided  with  the  exception 
of  the  uncorrected  line  length. 

Tape  D,  which  was  processed  some  months  before  tape  0,  has  less  data  in  its 
ancillary  blocks.  In  particular,  the  time  and  day/time  are  zero,  even  for  band  4. 
On  this  tape  the  EOV  flag  is  -27636. 

It  seems  possible,  therefore,  that  the  corrected  data  tapes  have  gone 
through  a  number  of  versions,  each  more  complete  than  the  previous  one. 

(c)  The  new  corrected  tape  has  almost  complete  ancillary  blocks  for  all 
bands.  The  scene  does  not  start  with  Sensor  No.1.  The  XSCAN  values  represent 
the  distance  along  the  swath  in  metres;  0,  80,  160  etc,  for  rows  1,  2,  3  etc. 

3.5.2  Pixel  start  and  end  numbers 

Each  MSS  data  record  has  3600  bytes  available  for  radiometric  data.  As 
mentioned,  the  data  is  stored  in  bytes  l8l  to  3780  for  band  4,  and  bytes  3  to 
3602  for  bands  5i  6  and  7- 

Normally,  considerably  less  than  36OO  radiometric  values  are  provided:  the 
actual  number  is  quoted  as  being  between  3000  and  3450  bytes,  and  in  practice  is 
about  3244  on  the  raw  data  tape.  The  actual  number  of  values  depends  upon  the 
interval  between  the  start-of-scan  and  end-of-scan  signals  in  the  satellite,  and 
seems  to  vary  very  little  over  a  considerable  period  of  time. 

All  values  are  stored  as  binary,  and  bytes  which  do  not  contain  radiometric 
values  are  set  to  zero,  such  values  occurring  both  before  and  after  the  radio- 
metric  values.  Since  zero  is  a  valid  radiometric  value  it  is  not  always  possible 
to  distinguish  the  start  and  end  of  data,  by  inspecting  the  bytes.  It  is  there¬ 
fore  necessary  for  the  system  to  specify  where  radiometric  values  start  and  end. 

All  of  the  tapes  contain  pixel  start  and  end  numbers  in  the  ancillary  block 
for  band  4  of  every  row  of  the  scene.  In  some  cases  (eg  tape  B)  the  ancillary 
blocks  for  bands  5>  6  and  7  also  include  start  and  end  numbers,  and  these  are 
always  the  same  as  those  given  for  the  corresponding  band  4.  However,  the  inter¬ 
pretation  is  different,  in  some  cases  for  every  band,  as  described  next,  where: 


SP  =  start  pixel  number  (bytes) 

FP  =  finish  or  end  pixel  number  (bytes) 

(a)  Raw  data  tape  A 


Start  at 
SP 

SP-178 

SP-178 

SP-178 


End  at 

FP  +  178  +  4 
FP  +  2 
FP 

(FP-2) 


The  end  for  band  7  is  deduced.  Since  the  radiometric  values  at  the  end  of 
the  line  were  zero  it  was  not  possible  to  determine  the  end  with  certainty. 

For  the  first  row  of  the  scene  SP  =  411  and  FP  =  3473  so: 


Band 

Start 

End 

No.  of  pixels  | 

4 

411 

3655 

3245 

5 

233 

3475 

3243 

6 

233 

3473 

3241  i 

7 

233 

(3471 ) 

(3239) 

The  values  of  SP 

and  FP  were  examined  for 

the  whole  of  the  tape.  SP 

remained  at  the  value  4l1  throughout  the  entire 

scene,  whilst  FP  alternated 

between  the  values  3473  and  3474, 

the  change  occurring  after  irregular  numbers  of 

lines,  eg  at  lines  19, 

25,  79,  85,  91  etc. 

(b)  New  corrected  data  tape  B 

Band 

Start  at 

End  at 

4 

SP 

FP 

3 

SP-178-2 

FP-178-2 

6 

SP-178-4 

FP-178-4 

7 

SP-178-6 

FP- 178-6 

For  the  first  row  of  this  tape  SP  =  4^4, 

FP  =  3753,  hence: 

Band 

Start 

End 

No.  of  pixels 

4 

454 

3753 

3300 

5 

274 

3573 

3300 

6 

272 

3571 

3300 

7 

270 

3569 

3300 

00 

r* 

O 
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This  tape  gives  the  'uncorrected  line  length'  as  3240,  for  band  4  row  1. 
This  does  not  refer  to  the  number  of  pixels  in  the  line  (which  is  seen,  above,  to 
be  3300)  but  corresponds  approximately  to  the  number  of  pixels  in  the  line  of  raw 
data  (see  (a)  above). 

The  values  of  SP  and  FP  were  found  to  decrease  progressively  throughout 
this  tape.  For  example,  SP  =  454  and  FP  =  3753  (for  all  four  ancillary  blocks) 
for  the  first  row  of  the  scene.  Subsequently,  SP  and  FP  both  decreased  by  one, 
at  about  every  20th  row  of  the  scene,  ending  the  scene  with  values  SP  =  348  and 
FP  =  3647.  In  every  case  the  difference  FP-SP  was  the  same,  3299,  ie  a  constant 
line  length  of  3300  pixels. 

The  decreasing  SP  and  FP  values  imply  that  the  scene  as  a  whole  is  a  para¬ 
llelogram  rather  than  a  rectangle,  this  being  a  result  of  geometric  correction 
for  Earth  rotation  below  the  moving  satellite.  It  is  normal  practice  in  RAE  to 
extract  the  maximum  possible  rectangle  from  this  parallelogram. 

(c)  Old  corrected  data  tape  C 


Band 

Start  at 

End  at 

4 

SP 

FP  +  178  +  8 

5 

SP-178 

FP  +  6 

6 

SP-178 

FP  +  4 

7 

SP-178 

FP  +  2 

first  row. 

SP  =  321 

,  FP  =  3415 

so: 

Band 

Start 

End 

No.  of  pixels 

4 

321 

3601 

3281 

5 

143 

3421 

3279 

6 

143 

3419 

3277 

7 

143 

3417 

3275 

In  the  older  corrected  tape,  the  uncorrected  line  length  is  not  given.  The 
number  of  pixels  for  each  band  reflects  the  varying  number  in  the  raw  data.  How¬ 
ever,  in  the  new  corrected  tape  (B)  the  number  of  pixels  has  been  made  constant 
(3300)  for  each  band. 

It  is  thus  apparent  that  there  is  an  older  (C)  and  newer  (B)  type  of 
corrected  data  tape,  which  need  different  interpretation,  at  least  concerning  the 
start  and  end  pixels. 
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ADDED  PIXELS 


It  has  been  shown  that  the  raw  data  scene  has  about  3240  pixels  per  row, 
whereas  the  corrected  scene  has  3300.  The  reason  for  this  is  that  extra  pixels 
have  been  inserted,  in  a  systematic  manner,  to  'correct'  the  geometry  of  the 
scene. 

There  are  two  major  distortions  in  an  MSS  scene  which  affect  the  distribu¬ 
tion  of  pixels  along  a  line.  Firstly,  the  oscillating  mirror  in  the  Landsat 
satellite  does  not  move  at  a  constant  speed,  and  this  gives  rise  to  'non-linear 
mirror  velocity  distortion'.  Secondly,  because  of  the  5*9  degree  field  of  view 
on  either  side  of  the  vertical,  there  is  the  so-called  'panoramic  distortion*. 
These  two  distortions  have  a  cumulative  effect,  which  is  to  'stretch  out'  the 
edges  of  the  scene,  in  relation  to  the  central  regions. 

To  compensate  for  these  distortions,  the  corrected  tapes  have  pixels  added 
to  them.  Pixels  are  added  only  infrequently  at  the  edges  of  the  scene  (ie  near 
the  ends  of  the  lines)  but  increasingly  more  are  added  near  the  centre,  for 
example  at  a  rate  of  one  extra  pixel  per  28  pixels.  This  effectively  ’stretches' 
the  central  regions,  so  that  the  whole  scene  is  linearly  compensated. 

4.1  Comparison  of  raw  and  corrected  tapes 

Information  about  the  added  pixels  can  be  obtained  by  comparison  of  raw  and 
corrected  tapes.  Each  line  of  a  corrected  scene  should  correspond  to  the  same 
line  of  the  raw  scene,  but  with  the  radiometric  values  transformed  by  means  of 
the  look-up  tables  and  with  added  pixels. 

When  the  two  tapes  were  studied  it  became  apparent  that  the  scenes,  although 
nominally  similar,  were  not  identical,  and  it  was  eventually  determined  that 
line  1  of  the  corrected  tape  corresponded  to  line  262  of  the  raw  data  tape.  With 
this  knowledge  it  was  then  possible  to  list  the  pixel  values  for  the  correspond¬ 
ing  lines. 

The  values  of  the  two  lines  would  not  be  expected  to  correspond  directly, 
because  the  raw  values  are  converted,  by  means  of  a  look-up  table,  to  corrected 
values.  When  this  conversion  was  taken  into  account,  the  two  lines  could  be 
compared.  It  then  became  apparent  that  extra  pixels  were  in  fact  being  added,  by 
duplication  of  preceding  values,  at  intervals  throughout  the  line. 

Following  this  discovery,  a  program  was  written  to  detect  all  pixels  which 
had  the  same  value  as  the  preceding  one.  Clearly  this  would  detect  both  added 
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pixels  and  other,  genuine  duplicates.  However,  by  repeating  the  process  for 
several  lines,  it  was  possible  to  detect  the  locations  of  the  added  pixels. 

When  this  program  was  run  on  50  lines  of  the  corrected  tape  B  (10  was 
insufficient)  it  showed  very  clearly  where  the  pixels  were  added.  At  the  ends  of 
the  line,  one  pixel  was  added  for  about  every  100  existing  ones,  and  at  the 
centre  one  pixel  per  30  was  added,  (ie  30  pixels  became  31),  with  79  pixels  added 
in  all. 

When  more  extended  tests  were  run,  it  became  apparent  that  the  sequence  of 
added  pixels  was  changed  at  intervals  throughout  the  scene.  Further  examination 
indicated  that  this  change  occurred  after  every  82  lines.  It  is  not  known 
whether  this  interval  is  common  to  all  new  corrected  tapes,  or  just  to  this  par¬ 
ticular  tape;  verification  would  be  a  lengthy  process. 

4.2  Study  of  older  corrected  tape 

An  older  corrected  data  tape  (C)  was  then  studied. 

When  the  program  described  above  was  run,  no  'added  pixels'  were  found. 
Further  studies  indicated  that  pixels  were  being  added  by  averaging,  rather  than 
by  duplication  of  existing  ones,  ie  each  added  pixel  was  assigned  a  value  which 
was  the  average  of  the  pixels  on  either  side  of  it. 

It  seemed  that  the  best  approach  would  be  to  inspect  the  look-up  table 
appropriate  to  the  band  and  sensor,  and  then  examine  the  line  for  radiometric 
values  not  in  the  look-up  table.  This  would  not  detect  all  added  pixels,  as 
sometimes  the  average  would  be  an  existing  value,  but  if  repeated  on  several 
lines  it  should  reveal  added  pixels  only. 

Another  program  was  prepared  to  perform  this  task,  and  this  successfully 
detected  the  added  pixels.  As  with  the  tape  B,  tape  C  shows  that  pixels  are 
added  more  closely  in  the  centre  of  the  line,  and  less  frequently  at  the  ends, 
hut  the  added  pixels  are  of  course  averaged  rather  than  duplicated  values. 

Adding  ’averaged*  pixel  values  is  not  acceptable  when  computer  interpreta¬ 
tion  is  to  be  done.  For  example  (and  for  simplicity  the  example  will  be  confined 
to  one  band,  -  although  it  is  equally  applicable  to  all  4  bands)  consider  two 
types  of  feature  on  a  scene,  feature  1  having  radiometric  values  of  say  10-20  and 
feature  2  having  values  of  90-100.  If  a  pixel  is  inserted  between  two  such 
features,  an  average  value  of  say  50-60  will  be  provided,  which  is  unlike  either 
feature  1  or  feature  2,  but  may  be  consistent  with  some  other  type  of  feature. 
Thus  adding  averaged  pixel  values  is  a  form  of  data  corruption. 
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4.3  Comments  on  added  pixels 

Geometric  correction  of  any  kind,  on  a  raster-scan  image,  must  inevitably 
give  rise  to  some  amount  of  data  degradation,  so  for  some  purposes,  the  raw,  and 
hence  distorted,  scene  is  preferable.  For  many  users,  however,  including  those 
performing  a  computer  analysis,  geometric  correction  may  have  some  advantages. 

This  study  has  shown  that  the  corrected  tapes  use  added  pixels  to  obtain  a  geo¬ 
metrically  corrected  scene,  and  that  the  'old'  method  was  to  add  averaged  values, 
whilst  the  'new'  method  is  to  insert  duplicate  values.  Both  methods  have  dis¬ 
advantages,  but  it  seems  likely  that  the  new  method  is  less  undesirable. 

This  study  has  shown  that  it  is  in  practice  possible  to  detect  the  locations 
of  the  added  pixels  with  both  old  and  new  tapes.  In  the  case  of  the  new  tape  the 
change  after  every  82  lines  makes  this  detection  very  laborious.  It  is,  however, 
practical  to  use  a  computer  to  detect  and  convert  the  added  averaged  pixels  of 
the  old  tapes  to  added  duplicate  pixels,  and  such  a  conversion  may  often  by 
advantageous,  for  image  analysis  purposes. 

5  PROGRAMS  TO  PROCESS  TAPES 

The  purpose  of  this  work  was  to  obtain  information  about  Landsat  MSS  data 
tapes,  so  that  'raw'  information  could  be  extracted  from  the  extensive  stock  of 
corrected  tapes  held  at  RAE.  During  the  course  of  the  work,  the  'new'  corrected 
format  became  the  standard,  and  consequently  demanded  a  new  program  to  read  it 
correctly. 

Consequent  upon  this  study,  therefore,  two  new  programs  have  been  prepared 
and  these  are  described  in  sections  5.2  and  5*5- 

5.1  Disc  image  files 

When  it  is  intended  to  perform  any  computer  analysis  of  a  Landsat  MSS  scene, 
it  is  convenient  to  read  the  CCT  and  form  a  set  of  image  files,  held  in  disc 
store.  This  is  done  chiefly  for  speed  of  access,  as  any  portion  of  a  disc  file 
can  be  accessed  in  much  less  than  one  second,  as  distinct  from  a  magnetic  tape 
where  access  can  take  several  minutes. 

It  is  convenient  to  form  four  equal  sized  disc  files,  of  bands  4  to  7 
inclusive.  Each  file  consists  of  a  'header'  followed  by  the  video  data  held  as 
a  raster  scan,  stored  in  sequence,  and  all  files  are  'in  register',  that  is  to 
say  they  correspond  to  exactly  the  same  portion  of  the  Earth's  surface.  Failure 
to  achieve  registration  will  of  course  falsify  any  inter-band  analysis. 
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Due  to  the  mode  of  operation  of  the  Landsat  Multi-spectral  Scanner,  the 
pixels  in  any  line  are  displaced  by  two  from  band  to  band  (see  Ref  1).  The 
specifications  mentioned  in  Refs  2  and  3  do  not  make  clear  which  pixels  are  in 
register.  During  the  development  of  the  programs  described  below,  it  has  been 
necessary  to  examine  the  disc  images,  to  check  for  registration,  and  to  modify 
the  program  accordingly.  Registration  is  most  easily  checked  by  selecting  some 
image  feature  which  is  conspicuous  in  all  four  bands,  and  ensuring  that  its 
location  is  identical  on  all  four  image  files. 

5.2  Program  to  convert  corrected  to  raw  data 

Since  'correction'  includes  both  the  conversion  of  values,  using  look-up 
tables,  and  the  addition  of  pixels  for  geometric  correction,  two  kinds  of 
•uncorrection'  are  possible:  these  being  of  course  the  removal  of  added  pixels 
and  the  restoration  of  the  original  radiometric  values.  It  would  be  possible  to 
write  a  program  to  perform  one  or  both  of  these  operations  for  the  'old'  format, 
and  another  program  to  do  the  same  for  the  'new*  format. 

It  was  finally  decided  to  prepare  a  program  which  could  operate  on  either 
old  or  new  tapes.  This  program  examines  the  date  of  production  of  the  tape, 
which  is  included  in  the  Landsat  header,  and  decides  whether  it  is  old  or  new 
according  to  whether  the  production  date  is  before  or  after  1  August  1979- 
Having  decided  this,  the  appropriate  conversion  is  selected. 

For  many  purposes  it  is  desirable  to  retain  the  geometric  correction,  and 
for  new  tapes  the  added  (duplicate)  pixels  are  acceptable.  However  the  added 
(averaged)  pixels  of  the  old  tapes  are  not  considered  acceptable,  and  thus  these 
are  replaced  by  duplicate  values.  In  the  case  of  both  old  and  new  tapes  the 
radiometric  values  are  then  restored  to  their  raw  values  by  operation  of  the 
look-up  tables  in  reverse. 

The  complete  program  is  named  IM.UNCOR,  consistently  with  the  naming  of  the 
other  image  programs  in  Space  Department.  The  listing  of  this  program  is  attached 
in  Appendix  A. 

Program  to  read  old  and  new  corrected  tapes 

Before  this  work  began,  there  was  a  program,  named  IM.FUC.MSS,  which  was 
used  to  read  the  old  (then  standard)  MSS  tapes  produced  in  Fucino  and  to  form 
image  files  from  the  data. 

A  new  program  has  now  been  prepared,  named  IM.FUC.MSS .2  which  incorporates 
both  the  facilities  of  the  earlier  program,  and  in  addition  the  ability  to 
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process  the  new  tapes.  As  with  IM.UNCOR,  the  program  itself  decides  whether  the 
tape  is  an  old  or  new  one,  and  proceeds  accordingly. 

The  opportunity  has  been  taken  to  incorporate  a  number  of  features  in  this 
program  which  make  ic  considerably  more  convenient  to  use  than  the  original  one, 
namely : 

(a)  The  program  determines  whether  the  character  code  used  is  EBCDIC  or 
ASCII,  instead  of  asking  the  user.  The  information  is  obtained  from  one  of  the 
process  flag  bits  in  the  Landsat  header  record. 

(b)  The  program  only  asks  for  one  disc  file  name,  and  suffixes  it  with 
4,  5»  6  or  7  according  to  the  band  in  question.  Likewise  it  only  asks  for  one 
'text'  for  the  file  header,  and  attaches  a  suffix.  Previously  the  user  had  to 
provide  four  file  names  each  with  text. 

(c)  The  file  names  and  texts  were  previously  requested  after  the  first 
pass  through  the  tape  had  been  done,  ie  about  30  minutes  after  the  start  of  the 
program.  Now  the  questions  occur  at  the  start  of  the  run,  so  that  no  waiting  is 
involved. 

(d)  The  user  can  specify  which  bands  he  wishes  to  use,  whereas  previously 
he  had  to  have  all  four.  This  is  advantageous  if  the  space  available  on  disc  is 
limited. 

(e)  The  user  can  specify  what  portion  of  the  image  he  wishes  to  use  - 
eg  lines  1000  to  2000,  whereas  before  the  whole  image  was  input. 

A  listing  of  this  program  is  given  in  Appendix  B. 

6  CONCLUSIONS 

A  detailed  study  has  been  made  of  a  raw  Landsat  MSS  CCT  and  severed  correc¬ 
ted  tapes,  both  old  and  new.  This  study  reveals  the  differences  between  the 
various  types  of  tape,  and  indicates  how  the  conversions  from  raw  to  corrected 
form  have  been  made. 

From  the  knowledge  gained  of  the  detailed  structure  of  these  tapes,  it  has 
been  possible  to  determine  how  to  convert  corrected  tapes  (of  which  RAE  Space 
Department  has  an  extensive  stock)  back  to  raw  data  form.  A  program  has  been 
written  to  do  this  for  both  old  and  new  corrected  tapes.  The  resulting  raw  data 
is  more  suitable  than  corrected  data  for  certain  computer  analyses. 

As  a  result  of  this  work,  a  new  version  of  the  tape  to  disc  file  input 
program  has  also  been  prepared.  This  has  several  operating  advantages  over  the 


existing  program  in  ease  and  versatility  of  use,  in  addition  to  being  capable  of 
processing  both  old  and  new  tapes,  whereas  the  existing  version  could  only  handle 
the  old  tapes. 
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LISTING  OF  PROGRAM  IM.UNCOR 


PACE  Bill 


C  PROGRAM  IH.UHCORi  TO  UNCORRECT  AM  OLD  OR  MEU  CORRECTED  FUCIMO  TAPE 


C  PROGRAM  IH.UNCOR'  TO  UNCORRCCT  AN  OLD  OR  MEU  CORRECTED  FUCIMO  TAPE 
C  FOR  OLD  TAPES  IT  CONVERTS  INTERPOLATED  ADDED  PIXELS  TO  DUPLICATES 
C  TNEN  (OLD  OR  NEW)  IT  CONVERTS  PIX  VALUES  BACK.  USING  INVERSE  LOOK-UP  TABLES. 
C  THE  OUTPUT  IS  ONE  TO  FOUR  IMAGE  FILES 

PARAMETER  H-37S8  /'•NO.  OF  BYTE8  PER  IMAGE  BLOCK 

LOGICAL  EOT /EOF. ERR. 0/ ASCII /NEW. FLAG /PRINTH 
INTEGER**  DNUM 

DIMENSION  I  BUFF  (  N>.ISTORE(  H)»MAHE<  IS  >.  ITEXTdi  ) 

DIMENSION  NBN< 4  >.NOB( 4. 4. 44  ).MBAK( 4. 4. 294 ) 

•  INSERT  SYSCOH >ASKE YS 

DATA  SUFFIX  /'4H434?/' 

DATA  NLINES/'SB/'  /'•NO  OF  LINES  TO  CHECK  FOR  ADDED  PIX 
DATA  IS.  IF/94.99/'  /'•USED  FOR  CALC  OF  1SP.IFP 
C 
C 
C 

URITEd.  IBBB) 

I IBB  FORHATt ' IN . UNCOR  IB-JAN-OB' > 

C 

21  Q-RNUN«A(' INPUT  TAPE  UNIT  NUMBER* .22. A*DEC. DNUN> 

IF  ( . NOT .0 >  GOTO  21 
I MUHI T“I NT 8< DNUM  ) 

IF  (IHUNIT.LT. I  OR.  1MUNIT . GT . 3)  GOTO  2B 
CALL  REU( INUNIT) 

C 

DO  3B  HBAHD-1.4 

3B  HBN(HBAND)-0  V*HO  BAND  IN  USE.. 

0-Y8H0SA('ALL  SANDS' . 9. ASNDEF > 

IF  ( . NOT  .0  >  GOTO  33 
DO  32  NBAHD-1/4 

32  HBH(HBAND)-! 

GOTO  41 

33  0-RNUHSA('N0U  MANY  BANDS' . 14.ASDEC. DNUM) 

IF  ( . NOT . 0 )  GOTO  33 
NBAND8-INT8(0NUM) 

IF  (NBANDS.LT. I  .OR. NBANDS .GT. 3  )  GOTO  33 
DO  39  J»l. NBANDS 

39  0«RNUNSA( ' BAND  NO' . 7 . ASOEC.DNUH  ) 

IF  ( .NOT.O)  GOTO  39 
<IB*INT8(  DNUM) 

IF  (JB.GE.4  .AND.  JS.LE.7)  GOTO  30 
37  URITEd.  1BB2) 

1BB2  FORMAT! ' REPEAT  THAT  NUMBER') 

GOTO  39 
30  JC-4B-3 

IF  (NBN( JC>. £9. I )  GOTO  37 
NBN(4C)>1 
39  CONTINUE 
C 

4B  LSTART-I  /'•START  FOR  ENTIRE  IMAGE 

NLINE-22BS  /'•NO.  OF  LINES  IN  ENTIRE  IMAGE 
LPI-NLINE  /'•END  LINE  FOR  ENTIRE  IMAGE 
0-YSN0SA(' ENTIRE  IMAGE' . 12. ASNDEF  ) 

IF  (0)  GOTO  48 

42  flaRNUMSA( *  STARTING  LINE' . 13. ASOEC.DNUH > 

IF  (.NOT.O)  GOTO  42 
L8TART-INTSCDNUH) 

IF  (L8TART.LT. I  .OR. L8TART .GT. 2284)  GOTO  42 
44  G*RHUHtA( ' NOU  MANY  LINES' . 14. ASOEC.DNUH) 

IF  (.NOT.O)  GOTO  44 
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C  PROGRAM  IN.UNCORI  TO  UNCORRECT  AN  OLD  OR  NEU  CORRECTED  FUCIHO  TAPE 


NL!NE-!NT8<DNUM> 

IF  INLINE. LT . 0  .OR.  HLINE . GT . 2286 >  COTO  44 
LPI-L START ♦ ML IHE-1  /‘LAST  LINE  OF  OUTPUT  I  MACE 
IF  (LPI.CT.220S)  GOTO  42 

READ  TEXT  AND  FILENAME  FROM  OPERATOR 
40  WRITE! 1. 1IR4> 

1114  FORMAT! 'ENTER  TEXT  MESSACE  FOR  IMAGE') 

READ! I. 1BB3 >ITEXT 
IIRS  F0RMATI4BA2 ) 

47  O-RNAMSA!' OUTPUT  FILE  NAME' . 10 , A9FUPP. NAME, 32 > 

IF  < . NOT .0 )  GOTO  47 

C 

C  READ  THROUGH  HEADERS,  HITH  OPTIONAL  PRINTOUT 
PRINTH-YSNOSA!' PRINT  HEADERS' . 13.PSNDEF) 

CALL  HTRE< INUNIT, IBUFF. 1S3B, EOT. EOF, ERR, NU)  /*READ  SLOCK  1 
IF  (PRINTN )  CALL  J8CPT! IBUFF)  /‘WRITE  JSC  HDR  TO  VDU 
CALL  NTRE! INUNIT. IBUFF, 193B, EOT. EOF, ERR, NU)  /‘READ  EOF 
IF  (  NOT. EOF)  WRITE! 1 , 1BBB  ) 

I  BBS  FORMAT! ' HOT  EOF  FOLLOUIHC  JSC  HEADER') 

CALL  NTRE< IHUNIT. IBUFF, 72B, EOT, EOF. ERR, NU)  /‘LANDS AT  HDR 
CALL  LAHPT! IBUFF. ASCII. NEU. PRINTH)  /‘LANDS AT  HEADER 
IF  (NEW)  WRITE! 1. IBIl  ) 

iaai  FORMAT! ' THIS  IS  A  ■NEU"  TAPE  -  PRODUCED  AFTER  JULY  79') 

CALL  NTREIINUN1T, IBUFF, 30B, EOT. EOF. ERR, NU)  /‘TRANSFORMATION  RECORD 
IF  IPRINTN )  CALL  TRANPT! IBUFF. ASC 1 1  )  /‘WRITE  TRANSFORMATION  RECORD 
DO  49  J-t.3  /‘DATA  RECORDS  4  TO  8 

CALL  NTRE!IHUNIT, IBUFF, 81B, EOT, EOF. ERR. HU)  /‘LOOK-UP  RECORDS 

IF  IJ.E8.3  .OR.  NBN! J  ) . NE . I >  GOTO  49 

JB-J‘3 

IF  !PRINTH )  CALL  LOOKPT! IBUFF. JB, ASCII  )  /‘WRITE  A  TABLE 
IF  !  NOT. PRINTH  .AND.  . NOT  ASCI  I >  CALL  EBCASC! IBUFF! 1 >. 81B.B  ) 

DECODE! 1336, till, 1BUFF)13T0RE  /‘TEMPORARILY  INTO  ISTORE 
Itll  F0RHAT!384! 14 )) 

DO  48  K-1.6 
KA-!K-1)‘64 
DO  482  L‘l , 236 

482  HBAK! J.K.L  )■-!  /‘SET  TO  -1  IF  HO  VALUE  EXISTS 
INO-8 

DO  484  L-1,84 
I N‘ IS TORE! KA‘L> 

H08! J , K, L  )■  IN 

IF-  ! INO. LT . 233)  NBAK! J. K, IN‘l )‘L- 1  /‘LOOK-UP  VALUES  FOR  UN-CORRECTINC 
I NO-IN 

4S4  CONTINUE 

48  CONTINUE 

49  CONTINUE 

CALL  MTRE! INUHIT. IBUFF. 133B. EOT, EOF, ERR, NU)  /‘READ  EOF 
IF  !. NOT. EOF)  WRITE! 1 . 1BB3  ) 

1BB3  FORMAT! ' NOT  EOF  AFTER  LOOK-UP  TABLES') 

C 

NWB-189B 

HUBl‘NUB‘t 

ISEHS-B  /‘START  AT  SENSOR  1  FOR  OLD  TAPES 
LINET-8  /‘COUNT  TEST  LINES 
DO  32  I-l.N 
32  ISTORE!! >‘B 

DO  IBB  LIME-1 .LP! 

ISENS-ISENS‘1 

IF  ! ISENS. GT . 6)  18EHS-1 
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C  PROGRAM  IH.UMCORI  TO  UNCORRECT  AN  OLD  OR  NEN  CORRECTED  FUCtHO  TAPE 


DO  OB  1BAND-1.4 

CALL  NTRE( I NUNI T . IBUFF. NUB , EOT . EOF, ERR , NW > 

IF  <  EOF  >  GOTO  998 
IF  (IBAND.NE.l)  GOTO  «B 
C  CHECK  ON  8ENSOR  NUMBER 

IS8-R8(IBUFF(4B>.8>  /*BYTE  117  IN  AHCIL  BLOCK 
IF  (NOT. NEN  .OR.  LINET . NE . I >  GOTO  93 
1888-188  /'•HEM  i  STARTING  SENSOR  -  SAVE  FOR  2ND  PASS 
ISENS-1SS 

S3  IF  (ISS.NE. ISENS)  GOTO  911 

IF  (LINE  HE  LSTART  .AND.  LINE.NE.LPI)  GOTO  SB 
ISP*2S4«R8( 1BUFF( 1 8) ,  8 )*RT( tBUFF( IS). 8) 

1FP-234«RS( IBUFF( IF).8)+RT( I8UFF( IF), 8) 

IF  (LINE. NE. LSTART)  GOTO  S4 

IST-ISP 

IFIN-IFP 

34  IF  (LINE.HE.LPI)  GOTO  SB 
I  ST 2" ISP 
IFIH2-IFP 

(B  IF  ( NEN >  GOTO  BB 

IF  (LINE. LT. LSTART)  GOTO  SB 

IF  (NBN( IBAHD). HE. 1)  GOTO  BB  /*USE  ONLV  SPECIFIED  BANDS 
IF  (LINET.  CE.NLINE8)  GOTO  BB  /'•CHECK  NLINE8  ONLY 
LINET-LINET+1 

DO  44  K-l.NMB  /'•UNPACK  BYTES  TO  MORDS . . . 

KA-NU81-K  /• . . .BACKHANDS# .  .  . 

KB-KA+KA  /•...USING  TNE  SANE  BUFFER 

IH-IBUFF(KA) 

IBUFF(K8-1)-R8(1N.B) 

44  IBUFF(KB)-RT(IH.8> 

Jl-ISP 

IF  (IBAND.NE.l)  Jl-Jl-178 

C  •OLD"  TAPES!  ISTORE-1  FOR  PIN  VALUES  NOT  IN  LOOK-UP  TABLE 
ITEST-IFP-ISPM81  /*NO  OF  PIN  IN  BAND  7  LINE 
KA-Jl*3-2*( IBAND-1  ) 

DO  78  K-1.ITE8T 

IPIX-IBUFF(K»KA)  /•8TART  AT  7TN  9TH  3RD  OR  1ST  PIN  - 
DO  48  L-1,44 

IF  (NOS( IB AND, I SENS/ L).EQ.IPIN)  GOTO  7B 
48  CONTINUE 
ISTORE(K  )*1 
7B  CONTINUE 
BB  CONTINUE 
IBB  CONTINUE 
C 

CALL  REN( INUNIT) 

IF  (NEN)  GOTO  128 
MR! TE( 1-1184) 

1184  FORNAT( ' ADDED  PIXELS  AT  POSITIONS!-') 

IOLD-B 

1TESTL-ITE8T+4  /•POSSIBLE  ISTORE  VALUES 
IF  (HEM)  !TESTL«1TEST 
DO  I1B  I “1 . ITESTL 
IF  ( I  STORE! I ) . EB .8  )  GOTO  118 
IDIFF-I-IOLD 
MRITE( 1. 1183)1- IDIFF 
UBS  F0RNAT(2I4> 

IOLD-I 

1 1 B  CONTINUE 

I-ITESTL-IOLD 


oo 
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C  PROGRAM  IH.UHCORi  TO  UNCORRECT  AN  OLO  OR  HEU  CORRECTED  FUCINO  TAPE 


Nft! T£< 1, 1 1 B7  > I 
111?  FORNAT(  tX< It) 

C 

121  IF  ( I ST2 . CT . I  ST )  IST-IST2 

IF  (IFIN2.LT  IFIN)  IF1N-IFIH2 
IF  (.NOT. HEW)  IEPL-IFIH-IST*181 
IF  (MEW)  IEPL-IF1M-IST-3 
WRI TE(  1.  UB6)NLINE<  IEPL 

lilt  FORMA T( *  IMAGE  SIZEi'.IS.'  ROWS  BY'. IS.'  COLS') 

C 

C  OPEN  OUTPUT  IMAGE  FILES  USING  TEXT  AND  NAME  OBTAINED  EARLIER  AND  IEPL 
ICHT£X-NLENSA(IT£XT.BB)*1  /*PDR  3186  P23-6 
ICHNAH-NLEHf A(NAHE/32  )♦ I  /'•FIRST  FREE  LOCATION 
CALL  1NITI 
DO  131  1-1.4 

IF  (NBN( I).NE .1 >  GOTO  13B 

CALL  HCHR*A(ITEXT.  ICNTEX.  SUFFIX,  I  >  /’•ADD  SUFFIX  TO  MESSAGE 
CALL  HCHR«A(NAHE/ICHNAM/8UFFIX. I >  /'•ADD  SUFFIX  TO  FILENAME 
CALL  OIWX( I.'ll'/MANE/IEPL.HLIHE. I  TEXT) 

131  CONTINUE 
C 

C  ••••••  NOW  FOR  SECOND  PASS  •••«•• 

PRINTH-. FALSE  . 

CALL  LHEADR(  1  NUN  IT  zPRINTH )  /'•NO  OUTPUT  TO  VDU  OH  SECOND  PASS 

C 

I SENS-8 

IF  (HEW)  I8EHS-ISSS-1 
DO  218  LIHE-l/LPI 
I  SENS- ISEHS+1 
IF  ( I 8ENS . CT . 6)  ISEH8-1 
DO  191  I BAND-1/  4 

CALL  NTRE( INUNI T . IBUFF. NWB. EOT . EOF. ERR. NN > 

IF  (EOF)  GOTO  991 

IF  (LINE .LT  .L8T ART  )  GOTO  191 

IF  (IBANO.EO.l)  18P-236*R3( IBUFFC I8).8)+RT( IBUFF ( IS  )/ 8) 

IF  (NBN( IB AND  ). NE . 1 )  GOTO  191 

DO  164  K-l.NWB  /'•UNPACK  BYTES  TO  WORDS... 

KA-NWB1-K  BACKWARDS/.  .  . 

KB-KA+KA  /'•...USING  THE  SAME  BUFFER 

IW-IBUFF(KA) 

IBUFF( KB-l >-R8( IW.8) 

164  IBUFF(K8 >-RT( IW/8) 

C  TRANSFER  THE  LINE.  ALTERING  ADDED  PIX  AND  UNCORRECTING 
K8HIFT-IST-1  /'•FOR  IBUFF  POINTER 
IF  (HAND.  NE.  1 )  KSHIFT-KSH1FT-178 
KIS-IST-18P-<6-2*( I8AND-1 )  ) 

ITEHP— 1 
FLAG-. FALSE. 

DO  l?l  K-l.IEPL 
IF  (NEW)  GOTO  166 
IF  ( ISTORE( K+KIS ). NE . 1 )  GOTO  166 
IF  (1TENP.HE.-1)  GOTO  169 
FLAG-. TRUE. 

GOTO  I7B 

166  ITENP-!BUFF(K*KSHIFT)tl 

ITENP-NBAK(  1BAMD.I8EN8.  ITEMP)  /'•UN-CORRECT  THE  PIX 
IF  (ITEHP. EG. -1  )  GOTO  971  /'•VALUE  NOT  IN  LOOK-UP  TABLE 
160  IF  (.NOT. FLAG)  GOTO  169 
IBIIFF(  K-l  )-ITEMP 
FLAG-. FALSE. 
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C  PPOCR AH  IH.UNCORi  TO  UNCORRECT  AN  OLD  OR  NEW  CORRECTED  FUCINO  TAPE 


1S9  I  BUFF (  K  >-I TEHP 
178  CONTINUE 

CALL  NLI2<  I  BAND/  18UFF  )  /‘•WRITE  OUTPUT  LINE  TO  FILE 
191  CONTINUE 
288  CONTINUE 
C 

CALL  RENT  INUNIT ) 

DO  258  1-1.4 

IF  <N8N< IJ.NE.l >  COTO  258 
CALL  CLOSEK I > 

258  CONTINUE 

STOP  i  1234  54 
C 

988  WRITEt  1.  1849 > ISS . 1  SENS 
1849  FORMAT! ' FAILS  SENSOR  CHECK.  ISS.  1SEHS s ' , 21 6 > 

COTO  995 
C 

978  IB-IBAND+3 

IF  (NOT. NEW)  WRITE! 1.187B1LINE, IB, K.IPIX 
1878  FORMAT! ' LINE' .15.' .  BAND' . 12. ' .  PIXEL  NO'. 15.'.  VALUE'. IS. 
•  '  NOT  IN  TABLE*  > 

C  COTO  995 

GOTO  1S8  /^TEMPORARY 
C 

988  URITE( 1. 1888) 

1888  FORHAT<* ERROR  988') 

COTO  995 
C 

998  NRITE( 1. 1898) 

1898  FORMAT! ' DOUBLE  FILE  NARK  REACHED') 

995  CALL  REH(IHUNIT) 

CALL  CALLI 
STOP  17777 
END 


SUBROUTINE  OIUX<  FUNI T . ET. NAHE, EPL . LP I , TEXT ) 
INTEGER  FUNIT.ET.EPL.LPI.TEXT(4R> 

DIHEHSION  NAHE(IS) 

INTEGER  EXTRA( 44  > 

INTEGER  STATE(1S).8TART(1S).CET(I6).CEPL( IS) 
INTEGER  CLPK IS).LNUH(IS) 

COHHON/IHCOH/STATE. START. CET . CEPL . CLP I , LNUN 


CHECK  UNIT  IS  CL08ED 

IF(STATE(FUNIT).EQ-8)  GOTO  4 
URI TEC  1 . 488 )  FUHIT 
488  FORHAT( ' •••  ERROR  •••'/ 

•  'ATTEMPT  TO  OPEN  1HACE  FOR  WRITING  ON  UNIT 

•  'WHEN  UNIT  IN  USE') 

CALL  CALLI 

CALL  EXIT 
4  CONTINUE 

C  CHECK  EPL  .CE.  8  AND  THAT  LPI  CE.  8 

I F<  EPL  .  CE . 1  .AND.  LPI  . CE .  1)  GOTO  28 
WRITE! 1. IBB >F UNIT 
188  FORHAT< ' •••  ERROR  •••'/ 

•  ' ATTEHPT  TO  OPEN  1HACE  FOR  WRITING  ON  UNIT 

•  'WITH  EPL  OR  LPI  LE88  THAN  ONE') 

CALL  CALLI 


'  .  12/ 


'  .  12/ 
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C  PROGRAM  IH.UNCORI  TO  UNCORRECT  AN  OLD  OR  HEM  CORRECTED  FUCINO  TAPE 


CALL  EXIT 
21  CONTINUE 

C  OPEN  FILE  AND  SET  EXTRA  INFO 
3  CALL  ODAHMX! FUNIT.  NAME  > 

START(FUNIT)  -  44 
CET(FUNIT)  -  ET 
CEPL(FUNIT)  -  EPL 
CLPK  FUNIT  >  -  LPI 
EXTRA< 1)  -  8TART<FUNIT> 

EXTRA<  2)  -  CET(FUNIT) 

EXTRA<  3  >  -  CEPL< FUNIT) 

EXTRA!  4)  ■  CLPK  FUNIT  ) 

C 

DO  1  1*1 <4B 

EXTRA(I  ♦  4>  -  TEXT(I) 

I  '  CONTINUE 

CALL  POSN<  FUNIT <  RBBBRI8  ) 

CALL  MTDAR<FUNIT<EXTRA< 44) 

LNUN< FUNIT)  -  1 

8TATE<  FUNIT )  -2  /•  OPEN  FOR  HRITING 

RETURN 

END 

C* •••••••••••••••••••••*•••••••• ••••*•*•*••••• 

SUBROUTINE  ODANMX< FUNIT. NAME > 

C  OPENS  A  DAN  FILE  FOR  UNITING 
•INSERT  SVSC0H>A*KEY8 
INTEGER  FUNIT 
LOGICAL  OK 
INTEGER  NAHEdS) 


OK  -  OPENS A< ASMRIT  +  ASOANF. NAME. 32 . FUNI T > 

I F< .NOT. OK)  CALL  EXIT 
CALL  P08N< FUNIT,  ■■■■III) 

RETURN 
END 

. . . 

SUBROUTINE  LHEADR< IUN1T .PRINTH  > 

C  READS  HEADER  BL0CK8 .  IF  PRINTH.  OUTPUTS  DETAILS  AND  ASKS  QUESTIONS 
INTEGER  BUFF< 1331)  /•LARGE  ENOUGH  FOR  JSC  HEADER 
LOGICAL  EOT <EOF<ERR<ASC II <NEN. PRINTH 
C 
C 

CALL  HTRE< I UNIT <  BUFF  < 193B<  EOT. EOF  <ERR<  NM  )  /*READ  BLOCK  1 
IF  < PRINTH )  CALL  JSCPT(BUFF)  /*MRI TE  JSC  HDR  TO  VDU 
CALL  NTRE< I  UN IT  < BUFF .  IS3B.  EOT.  EOF .ERR.  NM )  /*RE AD  EOF 
IF  (.NOT. EOF)  MR1TE< I <1112  ) 

1IB2  FORMAT! ' NOT  EOF  FOLLOMING  JSC  HEADER'  ) 

CALL  NTRE< I UN IT  <  BUFF. 721. EOT .EOF. ERR. NM)  /•LAHD8AT  HDR 
CALL  LANPT<  BUFF . ASCI 1 < HEM. PRINTH  >  /•LANDSAT  HEADER 
CALL  NTREdUHIT.BUFF.3SI.E0T.E0F.ERR.HM)  /•TRANSFORMATION  RECORD 
IF  (PRINTH)  CALL  TRANPT(BUFF< ASCI I )  /•MRITE  TRANSFORMATION  RECORD 
DO  3B  J-4.B  /•DATA  RECORDS  4  TO  8 

CALL  NTRE( 1UNIT. BUFF  <BII<EOT<EOF<  ERR. KM)  /*LOOK-UP  RECORDS 
IF  (PRINTH. AND. J. ME. 0)  CALL  LOOKPT( BUFF. J < ASC 1 1 )  /•MRITE  A  TABLE 
3B  CONTINUE 

CALL  HTRE< I UN IT . DUFF . 193B. EOT. EOF < ERR. NM )  /•READ  EOF 
IF  (.NOT. EOF)  MR!TE< 1 < 1BB3  ) 

1BB3  FORMAT! 'lIOT  EOF  AFTER  LOOK-UP  TABLES') 

RETURN 

o 
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C  PROGRAM  IN.UNCORi  TO  UNCORRECT  AH  OLD  OR  HEU  CORRECTED  FUCIHO  TAPE 


C* ••••••••••••••*•••••••••*•••••••••••••••*• ••*••*••*•••••••••••*•••**•• 

C  SUBROUTINES  TO  PRINT  HEADER  RECORDS  FROM  FUCIHO  LAHDSAT  TAPES 
C 

SUBROUTINE  JSCPT(IPT) 

C  SUBROUTINE  TO  PRINT  THE  J.S.C.  HEADER  < 3B6B  BYTES).  ALNAYS  EBCDIC  CHARS. 

DIMENSION  IPTC1  ).IBC 17) 

C 

C 

IB< 1>« LS<RT<IPT< 1369). 8)<B)  /*1ST  BYTE 

DO  18  1-1.8  /*NEXT  It  BYTES 
IC1-2*I 
IC2-IC1+1 

I B< IC1  )-LT< IPTC  I  +  13S9  >.  8> 

II  I B( IC2  )-LS( RTC IPTC !♦ 1 389 ) . 8 > . 8 > 

C 

CALL  EBCASCdB.  IS.  1)  /'•CONVERT  TO  ASCII.  UPPER  CASE 
CALL  EBCASC<IPT( 1444). 8.1 > 

HRITE< 1. 1BBB)  < I B( I > . I- 1 . 1 S > . I  PTC  1443 > . ( I  PTC  I > . I -1 444 , 1 447 . 

1BBI  FORNATC///' •••••••••  JSC  HEADER  •••••*••*• '// 

1  'SUN  ELEVATION  i'.8Al.'  MI LL I  RADS ' // 

2  'SUN  AZIMUTH  I'.BAl.'  MILL  I  RADS ' // 

3  'EARTH  ROTATION  l'.I8.'  MILLIRADS' // 

4  'SATELLITE  ALT  i'.4A2.‘  METRES'//) 

RETURN 

END 

C 

£••••••••••••••••••••••••••••*••••••••••••••••••••••••••••••••••••-••••• 

SUBROUTINE  LAHPTC IPT . ASCI  I . NEH. PR IHTH ) 

C  DETERMINE  ASCII.  HEU.  IF  REQD.  PRINT  LAHDSAT  HEADER  C 1441  BYTES) 

LOGICAL  ASCII. NEH. PRINTH 
DIMENSION  IPTC 1  ) 

C 

C 

ASCII-  TRUE . 

I- RTC 1PTCSB3).B)  /*FLAG  FOR  ASCII/EBCDIC 

IF  CI.E0.24B)  ASCII-. FALSE  .  /*24I-3SB  OCTAL  -EBCDIC  ZERO 

IF  C . NOT . A8CI I )  CALL  EBCASCC IPTC 1 >. 72B.B  )  /•CONVERT  TO  ASCII 
M-TPTC  484) 

M*  *RSC  M, 8)  /*F IR8T  BYTE 

H  1-.RT C  Ml . 4 )  /*KEEPS  4  BITS  ONLY 
N2-RTCM.4)  /•SECOND  BYTE.  4  CHARS  ONLY 
H-1I*H1+M2  /•ASSEMBLE  AS  A  NUMBER 
IYR-IPTC  489  ) 

I I- RSC IYR.8) 

1 1- RTC 11.4) 

12- RTC IYR. 4  ) 

I YR-1B«I 1*12 
NEH*. FALSE. 

IF  CIYR.CT.79)  HEU- . TRUE . 

IF  CIYR.E0.79  .AND.  M.CE.S)  NEH- . TRUE . 

IF  CPRINTH )  HRITEC 1. IBBfl)  C IPTC  I  >. I ■ 1 , 721  ) 

1111  FORMATC /'••••••••  LAHDSAT  HEADER  •••••••• ‘ ./ 18C 4BA2/ )/ > 

RETURN 

END 

. . . . . 

SUBROUTINE  TRANPTC IPT. ASCII)  /•PRINT  TRANSFORMATION  RECORD 
LOGICAL  ASCII 
DIMENSION  IPTC l > 

C 
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C  PROGRAM  I N . UHCOR i  TO  UNCORRECT  AH  OLD  OR  NEU  CORRECTED  FUCINO  TAPE 


C 

IF  (.HOT. ASCII)  CALL  EBCASC< IPTC1  >. 36B.B  ) 

URI  TEC  1.  1IIIX  IPT<  I  >»  1-1,368 > 

I IBB  F0RHAT( /'•*••••••  TRANSFORMATION  RECORD  ••••••*••' /9< /4BA2 >) 

RETURN 

END 

^••••••••••••••••••••••••••••••••••••••••••* . . 

SUBROUTINE  LOOKPT( IPT.HT.A8CII >  /'•PRINTS  RADIOMETRIC  LOOK-UP  TABLE  NT 
LOGICAL  A8C1I 
DIMENSION  IPT< 1 > 

C 

C 

IF  (.NOT. ASCII)  CALL  EBCASC< IPT< 1 ). 768.B > 

H-768 

IF  (NT.EB.8)  H-296  /'•THERMAL  BAND  8  HAS  2  SENSORS 


URITEO.  1BBR>HT.(IPT<  J).J-1.H) 

1111  F0RHAT<//'  RADIOMETRIC  LOOK-UP  TABLE  ••••••  BAND'. 12. 

•  6<  /4<  /32A2  )>  > 

RETURN 

END 

C****** ••••••••••••••••••••••••••••••••••••• •••••••••••••••••••••••••••• 


SUBROUTINE  EBCA8C( 1HORD.N. IUP)  /•CONVERTS  N  CHARS.  EBCDIC  TO  ASCII 
C  STARTING  FROM  IUORD 
C  ) 

INTEGER  EBTOAS. EBCNAR.ASCHAR 
INTEGER  IUORO<  1  ).LCNAR(  2).  <1.  M.  I  CHAR.  L 
INTEGER  E8T0AS( 27) 

DATA  EBT0AS/HR6619. H4B24B. >899241.  >177642. i B79643 , > 149244. 
l  >866249. >898244. >176647. >844698. IB966SI. I B96292 . I B47293 . I I696S4. 
i  II6B299. 1149696. H6I697. H79272. >896673. > B46274 . > B7727S .» 167276. 
l  >867477. > 176381. > 166737. >122612. >899644/ 

C 

C 

IF  (N  .LE.  DCOTO  191 
DO  18  J-l.N 

SPLIT  IUORD(J)  INTO  TOO  SINGLE  CHARACTERS  LCHAR(I)  4  LCHARC2) 

LCHAR(l)-  R8( (AHD<  IH0RD</l).  >  177488)  ).B) 

LCHAR(  2 )■  AND<IHORD<<I  ).  1377) 

DO  28  M-1.2 
ICHAR-LCMAR(M) 

L«1 


8EE  HNETNER  ICNAR  IS  ALPHANUMERIC 


IF  (ICNAR  .LT.  >361 
ICNAR-  ICNAR-ilM 
GOTO  120 

.OR. 

ICNAR  . 6T . 

>371 >COTO 

211 

2IR 

IF  (ICNAR  .LT.  >311 
GOTO  12B 

OR 

ICNAR  GT. 

>3U)G0T0 

211 

211 

IF  (ICNAR  .LT.  >321 
ICNAR-  ICNAR- >017 
GOTO  121 

.DR. 

ICNAR  .GT. 

>331 )COTO 

221 

221 

IF  (ICNAR  .LT.  >342 
ICNAR-  ICNAR- >117 
GOTO  121 

.OR. 

ICNAR  .GT. 

>391 )GOTO 

231 

231 

IF  (ICNAR  .LT.  >201 
ICNAR-  ICNAR*  >141 

.OR. 

ICNAR  GT. 

>211  )COTO 

241 

GOTO  13B 
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C  PROGRAM  IH.UNCORi  TO  UNCORRECT  AN  OLD  OR  NEU  CORRECTED  FUCIHO  TAPE 


248 

IF  (ICHAR 

.LT. 

1221 

OR. 

ICHAR 

GT. 

1231 >COTO 

2SB 

ICHAR-  I CHAR*  >131 

GOTO  13B 

2SB 

IF  (ICHAR 

•  LT. 

(242 

OR. 

ICHAR 

GT. 

>251 )GOTO 

26B 

I  CHAR"  I  CHAR* 1121 
GOTO  130 

SEE  HHETHER  ICHAR  IS  IN  THE  ARRAY  EBTOAS 

2SB  EBCHAR-  RS( ( AND( EBT0AS( L  ) «  >1774BB))>8) 

IF  < ICHAR  . EO.  EBCHAR)  GOTO  4B 
L"L*1 

IF  < L  LE.  27  )C0T0  2tfl 
GOTO  SB 

48  ASCHAR"  AN0< EBT0AS< L >, i 377  ) 

ICHAR  -  ASCHAR 
GOTO  12B 

SET  ICHAR  TO  NUL  CHARACTER 

SB  ICHAR  -  l 2BB 
GOTO  12B 

IF  IUP  -  1  UPCASE  ICHAR 

I3B  IF  (IUP  HE.  DGOTO  12B 
ICHAR"  ICHAR- iB4B 
12B  LCHAR(H)-  ICHAR 
2B  CONTINUE 

RECOMBINE  LCHAR< 1 )  AND  LCHARC 2 )  INTO  A  SINGLE  WORD  (IMORD(J)) 

I UORD( J  >■  0R(  L8( LCHARI 1  )»  8  ). LCH ARC  2  )  ) 

IB  CONTINUE 
1SB  RETURN 
END 
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LISTING  OF  PROGRAM  IM.FUC.MSS.2 


FhCE  BSBl 


C  PROGRAM  IH.FUC.MSS.2 


C  PROGRAM  IH.FUC.MSS.2 

C  TO  RE AO  A  FUCINO  MAGTAPE <  IN  OLD  OR  MEN  FORMAT.  TO  IMACE-FILES 
PARAMETER  M-3788  /*HO  OF  BYTES  PER  IMAGE  BLOCK 
PARAMETER  NBB-1898  /*NO.  OF  BOROS  PER  IMAGE  BLOCK 

LOGICAL  EOT . EOF  »  ERR.  0 . HEW . PRINTH 
INTEGER*^  DNUH 

DIMENSION  I BUFU!  NBB  ) • IBUFB( H >. I  STORE! H  >, I  TEXT! 4B >. HBN( 4  > 

DIMENSION  NAME! 1G> 

4INSERT  SYSCQM )A4K£YS 

DATA  SUFFIX  /4N45G7/' 

C 

C 

C 

BRITE! 1, 10B0) 

I BBB  FORMAT! ' IN . FUC. MSS .2  11-JAN-8B'> 

C 

2B  0-RHUH4A!' INPUT  TAPE  UNIT  NUMBER' . 22 . A*DEC. DNUH ) 

IF  <  NOT  Q  >  GOTO  2fl 
INUNI T”I NTSIDNUN ) 

IF  ! I MUN 1 T . LT  .8  OR.  INUHIT.GT  3>  COTO  28 
CALL  REB! IHUNIT > 

C 

00  38  J-1,4 

38  HBN! J  >»B  /'♦NO  BAND  IN  USE.. 

0-YSNC4A!' ALL  BANDS' / 9, ABNDEF > 

IF  ! . NOT .8  )  COTO  33 

DO  32  J»  1 . 4 

32  HBN! J )-l 
N BANOS-4 
GOTO  48 

33  0=»RHUH»A! 'HOB  MANY  BANDS' , 1 4 , Af DEC. DNUH > 

IF  ! . NOT  . Q  )  GOTO  33 
NBANDS-INTSIDNUH) 

IF  ! HBANOS  LT. 1  . OR . NBANDS . GT . 3 >  COTO  33 
DO  39  J-l, NBANDS 

35  0-RNUH4A! 'BAND  HO' . ? . A»DE C , DNUH  ) 

IF  <  NOT .8)  COTO  35 
JB-INTSIONUM) 

IF  IJB.GE.4  .AND.  J8.LE.7)  COTO  38 
37  BRITE! 1. 1BB2> 

1882  FORMAT! 'REPEAT  THAT  NUMBER' > 

GOTO  35 
3B  JC-JB-3 

IF  ! NBN! JC  >  .  EG . 1  )  GOTO  37 
NBN! JC  )-l 

39  CONTINUE 
C 

48  LSTART-l  /'♦START  FOR  ENTIRE  IHACE 

NLINE-2286  /♦NO.  OF  LINES  IN  ENT  IRE  IMAGE 
LPI-NLINE  /'♦END  LINE  FOR  ENTIRE  IMAGE 
0-YSH04A! 'ENTIRE  IMAGE' . 1 2 . A4NDEF  ) 

IF  !0>  GOTO  5B 

42  Q-RNUMtA! *  STARTING  LINE'/ 13/AfDEC/DNUM) 

IF  !  NOT .0)  GOTO  42 
LSTAR TAINTS! DNUH  > 

IF  !LSTART . LT . 1  . OR . L8T ART  CT . 2284 >  GOTO  42 
44  QaRNUNtA! ' HOB  MANY  LINES' . 14 .A»DEC, ONUM > 

IF  ! .NOT. 8)  COTO  44  c 

NLINE* INT8! ONUM  >  - 

IF  INLINE. LT. 8  OR.  NLI NE . GT . 2284 >  COTO  44  ° 

J 
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C  PROGRAM  IH.FUC.NSS.2 


IF  INLINE. EO.B)  GOTO  4B  /*FOR  HEADER  PRINTOUT  ONLY 
LPI*L START ♦NLINE-1  /'•LAST  LINE  OF  OUTPUT  INACE 
IF  !LP1.CT.2284)  GOTO  42 

READ  TEXT  AND  FILENAME  FROM  OPERATOR 
SI  WRITE! 1/ IBI4) 

1114  FORMAT! ' ENTER  TEXT  MESSAGE  FOR  IMAGE') 

READ!  1  .  IMS  >1  TEXT 
IMS  FORMAT! 4IA2  > 

32  fl-RHANSA!' OUTPUT  FILE  NAME '. 14. AtFUPP. NAME. 32 ) 

IF  ! . NOT .0 )  COTO  32 
C 

C  READ  THROUGH  HEADERS/  WITH  PRINTOUT 
40  PRINTH- .  TRUE .  /'•PRINT  THE  HEADER 

CALL  LNEADR! I MUN1T. PRINTH .NEW)  /•-!  FOR  PRINTOUT 
IF  INLINE.  EO.B)  STOP  H2343  /•HLIHE-B  FOR  HEADER  PRINT  ONLY 
C 

C  READ  TAPE  ONCE  TO  DETERMINE  MAXIMUM  START  AND  FINISH  COLUMNS 
DO  IBB  LINE-1. LPI 
DO  9B  IBAND-1.4 

CALL  NTRE! IHUN1 T . IBUFW. NWB .EOT . EOF. ERR.HW > 

IF  ! EOF )  COTO  99R 
IF  IIBANO.NE.l)  GOTO  9B 
IF  ILINE .NE .L START )  GOTO  PB 
IST-234*RS! ISUFW!34 ), 8>*RTT IBUFW! 34 >.8 > 

IFIH-234*RS! IBUFW! 33 )«8  >*RT! IBUFW! 93  >. 8> 

PR  IF  ILINE. NE.LPI)  COTO  9R 

IST2-234«RS! IBUFW! 34  >.8  >*RT! IBUFW!  34 >. B> 

I FIH2-234*RS! IBUFW! 33 ). 8)*RT! IBUFW! 33). 8 ) 

9B  CONTIHUE 
IBB  CONTINUE 

CALL  REW! IHUNIT ) 

IF  IIST2.CT.I8T)  IST-1ST2 
IF  IIFIN2.LT. IFIN)  IFIH-IFIN2 

IF  ! NEW >  IEPL-IFIN-IST-3  /'•NO.  OF  PIX  IN  LINE  FOR  OUTPUT  IMAGE 
IF  !. NOT. NEW)  IEPL-IFIN-ISTM81 
C 

C  OPEN  IMAGE  FILE8  FOR  WRITING  USING  TEXT  AND  NAME  OBTAINED  EARLIER 
CALL  INITI 

ICHTEX-NLEN4A! I  TEXT. 8I)  +  1  /*PDR  31B4  P23-4 
ICHNAM-NLEN4AINAME.32  )♦!  /•POSITION  FOR  SUFFIX 
DO  IBS  I BAND* 1. 4 
IF  INBN! IB AND  >. HE . 1 )  GOTO  IBS 

CALL  NCNRSAIITEXT. 1CNTEX. SUFFIX. IBAND)  /«PDR  31RS  P23-P 
CALL  NCNR8AINANE.ICNNAN. SUFFIX. IBAND)  /'•ADDS  SUFFIX 
CALL  OIWX! IBAND. Ml' .NAME. 1EPL.NLINE.ITEXT) 

IBS  CONTINUE 

WRITE! I. 1BB4 INLINE. 1EPL 

1BB4  FORMAT! 'IMAGE  SIZE  I '.IS.'  ROWS'. IS.'  COLUMNS') 

C 

PRINTH-. FALSE. 

CALL  LNEADR! IMUN1T.PR1HTN. NEW)  /*N0  OUTPUT  TO  VDU  ON  SECOND  PASS 
C 

DO  2BB  LINE-1. LPI 
DO  19B  IBAND-1.4 

CALL  NTRE! INUNIT. IBUFW. NWB. EOT. EOF. ERR. NN > 

IF  <EOF)  GOTO  990 

IF  <LINE .LT .LIT ART )  GOTO  I9B 

IF  (NBN! IBAND >.NE. 1)  GOTO  19R  /•IGNORE  THIS  BAND 
DO  HR  KM. NWB 
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C  PROGRAM  IM  .FUC  .M88. 2 


KK*K»K-1 

IBUFB* KK )-RS( IBUFU*K).8> 

111  IBUF8( KK*1 )"RT( IBUFU*  K)<8  ) 
K8NIFT-I8T-I 

IF  (IBAND.NE.l)  K8MIFT-KSHIFT-178 
00  13M  K-l. IEPL 
I  STORE* K  >* I BUFBC K»KSHIFT) 

131  CONTINUE 

CALL  ULI2* I  BAND. I8T0RE) 

191  CONTINUE 
21 ■  CONTINUE 
C 

CALL  REU* IHUMIT ) 

00  2SI  I>1.4 

IF  <NBN( I > . NE . 1 >  GOTO  29B 
CALL  CL08EKI) 

29B  CONTINUE 

STOP  H23498 
C 

99B  UNITE* 1.1B9B) 

1B9B  FORMAT* ' DOUBLE  FILE  MARK  REACHED') 
CALL  REH* IHUNIT ) 

STOP  >7777 
END 


SUBROUTINE  01 HX*  FUNIT  .ET« NAME. EPL .LPI. TEXT) 
INTEGER  FUNIT <ET<EPL<LP I. TEXT*  4B> 

D1NEHS10N  NAME* 18) 

INTEGER  EXTRAC44  )  . 

INTEGER  STATE* 18  )<  START* 18).CET*18).CEPL( 18) 
INTEGER  CLPI*18).LNUH*18) 

COHNON/IHCOM/STATE.START. CET . CEPL . CLPI . LHUH 


C  CHECK  UNIT  18  CLOSED 

IF(STATE(FUNIT).EQ.B>  GOTO  4 
WRITE* 1.4BB)  FUNIT 
4BB  FORMAT* '•••  ERROR  ••••/ 

•  'ATTEMPT  TO  OPEN  IMAGE  FOR  UR 1  TING  ON  UNIT  '.12/ 

•  'UHEH  UNIT  IN  USE') 

CALL  CALLI 

CALL  EXIT 
4  CONTINUE 

C  CHECK  EPL  .CE.  ■  AND  THAT  LPI  CE.  ■ 

IF*  EPL  GE. 1  .AND.  LPI  .GE.  1)  GOTO  21 
UNITE* 1. IBB  )FUNIT 
IBM  FORMAT*'***  ERROR  ••*'/ 

•  'ATTEMPT  TO  OPEN  IMAGE  FOR  UR1TING  ON  UNIT  '.12/ 

•  'UITN  EPL  OR  LPI  LESS  THAN  ONE') 

CALL  CALLI 

CALL  EXIT 
2B  CONTINUE 

C  OPEN  FILE  AND  SET  EXTRA  INFO 
3  CALL  OOANNX* FUN I T«  NAME ) 

START* FUNIT)  -  44 
CET(FUNIT)  -  ET 
CEPL* FUNIT)  *  EPL 
CLPI* FUNIT )  •  LPI 
EXTRA* 1)  •  START(FUNIT) 

EXTRA* 2)  •  CET* FUNIT) 
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C  PROGRAM  IN.FUC.NSS.2 


EXTRA( 3 >  -  CEPl( FUNIT  ) 

EXTRA!  4  >  -  CLPI(FUNIT) 

C 

00  1  I-1.4B 

EXTRACI  ♦  4)  ■  TEXT<  I  > 

1  CONTINUE 

call  posn< funit; aaaiiia) 

CALL  UT0AM(FUNIT;EXTRA;44) 

LNUM( FUNIT  >  >  1 

STATE! FUNIT  >  >2  /•  OPEN  FOR  URITINC 

RETURN 

END 

C* ••••**•••****••*•*•*••••••*•••••**• ••*•••••• 

SUBROUTINE  ODANWX! FUNIT, NAME > 

C  OPENS  A  DAM  FILE  FOR  WRITING 
SINSERT  SYSC0H>A4KEY8 
INTEGER  FUNIT 
LOGICAL  OK 
INTEGER  NAME! 18  ) 


OK  -  OPEN* A! A4WR1T  ♦  ASOANF, NAME, 32 , FUNI T > 
IF< .HOT. OK)  CALL  EXIT 
CALL  P08N< FUNIT,  BBBBBBB) 

RETURN 

END 


SUBROUTINE  LHEADR! IUN1T, PR1NTH, NEW) 

READS  HEADER  BLOCKS.  IF  PRIHTH  IS  TRUE,  PRINTS  DETAILS 
INTEGER  BUFF!  1S3B >  /'•LARGE  ENOUGH  FOR  JSC  HEADER 
LOGICAL  EOT, EQF. ERR, A8CII.HEU. PRIHTH 


CALL  MTRE! IUNIT, BUFF. IS3a, EOT. EOF, ERR, NH>  /*REAO  BLOCK  I 
IF  ! PRIHTH )  CALL  J8CPT!BUFF)  /•WRITE  JSC  HDR  TO  VDU 
CALL  MTRE! IUNIT, BUFF. 193R, EOT, EOF, ERR, NW >  /«READ  EOF 
IF  !. NOT. EOF)  WRITE!!, 1B82) 

IRB2  FORMAT! *  NOT  EOF  FOLLOWING  JSC  HEADER') 

CALL  MTRE! IUNIT, BUFF, 72B, EOT. EOF, ERR, NW)  /*LAHDSAT  HDR 
CALL  LANPT! BUFF , ASCI  I , NEW , PRIHTH)  /•DETERMINE  ASCII.  NEW 
CALL  MTRE! IUNIT, BUFF. 3SB, EOT, EOF. ERR, NW)  /•TRANSFORMATION  RECORD 
IF  !PR1NTN )  CALL  TRANFT(BUFF,A8CI1 >  /*WR1TE  TRANSFORMATION  RECORD 
DO  3R  J«4,8  /«DATA  RECORDS  4  TO  S 

CALL  MTRE! IUNIT. BUFF. SIR. EOT, EOF, ERR, NW)  /*LOOK-UP  RECORDS 
IF  (PRIHTH . AND. J . HE . 8  )  CALL  LOOKPT! BUFF, J , ASC 1 1 >  /•WRITE  A  TABLE 
3R  CONTINUE 

CALL  MTRE! IUNIT. BUFF. 1338, EOT, EOF. ERR, NW )  /*READ  EOF 
IF  (.NOT. EOF)  WRITE! I ,1813 ) 

IBB3  FORMAT! 'NOT  EOF  AFTER  LOOK-UP  TABLES') 

RETURN 

END 

»••••••••••••••••••••••••••*••*•»••••«••»••••••••**••*•••••••••*••••••• 

SUBROUTINE  JSCPT(IPT) 

SUBROUTINE  TO  PRINT  THE  J.8.C.  HEADER  ! 3RSB  BYTES).  ALWAYS  EBCDIC  CHARS. 
D1HENSI0N  1PT(1 ),I8! 17) 


I B( 1 >*LS(RT( IPT( I3SB  >,8 ), B  )  /«18T  BYTE 

DO  IB  1*1,8  /AMENT  14  BYTES 
ICI-2«I 


PACE  BUS 


C  PROS RAH  IN.FUC.H8S.2 


IC2-IC1+1 

IB<  1C1  )-LT<  IPT<  IM3t9).B) 

IE  I B( IC2  )*LS( RT< IPT< I* 1 3*9), 8), 8  ) 

C 

CALL  EBCASCCIB,  It,  1)  /•CONVERT  TO  ASCII ,  UPPER  CASE 
CALL  EBCASCIIPTI 1444), 8.1) 

NRI TE< 1, 1BIB >  <I8< I), 1*1. lt),IPT< 1443), ( I PT< I >.1-1444,144?) 
FORHAT< ///'•••••••*•  JSC  HEADER  *•*•*••••*'// 

1  'SUN  ELEVATION  I'.SAl.'  HILL  I  RADS' // 

2  'SUN  AZIMUTH  i'.BAI,'  HILL1 RADS' // 

3  'EARTH  ROTATION  I'.IS. '  MILLIARDS' // 

4  'SATELLITE  ALT  i',4A2.'  METRES'//) 

RETURN 

END 

C 

^••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••••t*»«»* 

SUBROUTINE  LANPT( IPT, ASCI  I ,MEH, PRINTH) 

C  DETERMINE  ASCII,  MEM .  IF  REBD.  PRINT  LANDSAT  HEADER  (144B  BYTES) 
LOCICAL  A8C1I, HEM, PRINTH 
DIMENSION  IPT< 1 ) 

DOUBLE  PRECISION  CHASC.CHEBC 
OATA  CHASC , CHEBC/8HA8C1 1  . BHEBCDIC 

C 
C 

I  *RT( IPT(S89),B)  /*FLAG  FOR  ASCII/EBCDIC 

IF  < I .NE .24B)  GOTO  21  /*IT'S  HOT  EBCDIC  (24B-3&B  OCTAL-EBCDIC  B 
ASCII*. FALSE. 

CALL  ERCA8C<IPT<  D.72B.B)  /•CONVERT  TO  ASCII 
IF  (PRINTH  )  MRITE(  1,  1BBDCHEBC 
IBB I  FORNAT<' MAGTAPE  CHARACTERS  ARE  ' , A8 ) 

GOTO  38 

21  ASCII* .TRUE . 

IF  (PRINTH)  MRITE< 1, 1BB1 )CMASC 
3B  M*IPT<  484  > 

NI-RS(H.S)  /-FIRST  BYTE 
N1*RT( Ml ,4  >  /»KEEPS  4  BITS  ONLY 
N2*RT( M. 4)  /•SECOND  BYTE.  4  CHAR8  ONLY 
N*1R*HI+M2  /•A88EM8LE  AS  A  NUMBER 
IYR*IPT<  483 ) 

I  t*RS<  IYR,  8  ) 

I t"RT< II, 4  > 

I2*RT< IYR, 4  ) 

I YR*1B*I 1-12 
MEM*. FALSE. 

IF  (IYR.GT.79)  HEM*. TRUE. 

IF  (IYR. EB . 79  .AND.  N.GE.B)  HEM*. TRUE. 

IF  (1ACT.LT.8)  MRITEt 1, 1RBR)  < IPT< I  ), 1*1 . 72B) 

1BBB  FORMAT< /'••••••••  LANDSAT  HEADER  ••••••••', /1B(4BA2/)/ ) 

RETURN 

END 

. . . 

SUBROUTINE  TRANPT( IFT, A8CI I >  /«PRINT  TRANSFORMATION  RECORD 
LOGICAL  ASCII 
0INENS10N  IFT<1> 

C 

C 

IF  (.NOT. ASCII)  CALL  EBCA8C< IPT( 1 >. 36B.S > 

MRITEI I.  IIIEXIFK  I ),  1*1, 3SI  > 

FORNATt /'••••••••  TRANSFORMATION  RECORD  •••••••••' /9( /4BA2 >> 

RETURN 
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C  PROGRAM  IH.FUC.HSB.2 


END 

C*  «*•••••*•••• . ••• . . . . 

SUBROUTINE  L00KPT( IPT. NT.  ASCII >  /'•PRINTS  RADIOMETRIC  LOOK-UP  TABLE  NT 
LOGICAL  A8CII 
DIMENSION  IPT< 1 ) 


IF  (.HOT. ASCII)  CALL  E8CASC(  IPT( 1 >, 818.8 > 

H-768 

IF  (NT.EO.S)  M-29S  /•THERMAL  BAND  8  HAS  2  SENSORS 
URI TE< 1. 1BBB>HT.(IPT( J>.J»1.N> 

IBBB  FORMAT!//'  RADIOMETRIC  LOOK-UP  TABLE  •*•»••  BAND'. 12. 
•  6(/4(/32A2 >>> 

RETURN 

END 


SUBROUTINE  E6CA8CI IBORD.H. IUP)  /-CONVERTS  N  CHARS.  EBCDIC  TO  ASCII 
STARTING  FROM  IUORD 

INTEGER  EBTOAS. EBCHAR < ASCHAR 
INTEGER  IHORD< 1 ).LCHAR(  2). J. H. ICNAR.L 
INTEGER  E8T0A8<  27) 

DATA  EBTOAS/>BB6619. > 848248. IB99241 . >877642. I B73643. I 849244. 

I  i B 64249 . i B3B246 . i 876647. i I4669B. >896631 . >896292. > 847293. *869694. 
i  H6B253. IB43636.IB6B637, H79272. IB34473. i 844274 . i B77279 .* B67276. 

I  1B67677. IB763BB.1B44737, >822612. >899644/ 


IF  (N  ,LE.  B)COTO  198 
DO  IB  J-l.N 

SPLIT  IMORD(J)  INTO  TOO  SINGLE  CHARACTERS  LCHAR(l)  6  LCHAR(2> 

LCHAR< 1 >■  R8( <  AND( IMORDO >.)  177411 >>. 8 > 

L CHART  2  >■  AND<  IHORD(  J  ).  >377) 

DO  28  M-1,2 
ICHAR-LCHAR(H) 

L-l 


8EE  MNETNER  ICHAR  IS  ALPHANUMERIC 


IF  ( ICNAR  .LT.  >368 
ICHAR-  ICHAR- >188 
GOTO  128 

.OR. 

ICHAR 

.GT. 

>371 >COTO 

288 

288 

IF  (ICHAR  .LT.  >381 
GOTO  128 

.OR. 

ICHAR 

.GT. 

•311 >COTO 

218 

218 

IF  (ICHAR  .LT.  >321 
ICHAR-  ICHAR- >887 
GOTO  128 

.OR. 

ICHAR 

■  GT. 

•  33D60T0 

228 

228 

IF  (ICHAR  .LT.  >342 
ICHAR-  ICHAR- >817 
GOTO  128 

.OR. 

ICHAR 

•  GT. 

>391 >60T0 

238 

238 

IF  (ICHAR  .LT.  >281 
ICHAR-  ICHAR- >148 
GOTO  138 

.OR. 

ICHAR 

.GT. 

>211 >COTO 

248 

248 

IF  (ICHAR  .LT.  >221 
ICHAR-  ICHAR- >131 
GOTO  138 

.OR. 

ICHAR 

■  GT. 

>231 JCOTO 

298 

298 

IF  (ICHAR  .LT.  >242 
ICHAR-  ICHAR- >121 

.08. 

ICHAR 

.GT. 

>291 >GOTO 

268 
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C  PROGRAM  IH.FUC.N88.2 


GOTO  131 

SEE  UHETHER  I CHAR  IS  IN  THE  ARRAY  EBT0A8 

RAM  EMC  MAR*  RS(<AMO(E*TOAS(L><  il7?4**)><8> 

IF  < ICHAR  . EO.  EICNAR)  GOTO  41 
L*L+1 

IF  (L  LE.  27>COTO  28* 

GOTO  SR  , 

4*  ASCHAR*  AN0<E0T0A8<D<  1 37?) 

I CHAR  -  ASCHAR 
GOTO  12* 

SET  1CHAR  TO  NUL  CHARACTER 

S*  ICHAR  *  1211 
GOTO  12* 

IF  1UP  -  I  UfCASE  ICHAR 

IX*  IF  UUP  .HE.  DGOTO  12* 

ICHAR*  I CHAR- i *4* 

12*  LCHAR<  N>*  ICHAR 
2*  CONTINUE 

RECONSIHE  L CHART  1 >  AMO  LCHAR<2)  INTO  A  SINGLE  UORO  <IN0R0<J>> 

IM0RD< J >■  OR<LS<LCHAR<I >.  S  >,  LCMAR<  2  >  > 

I*  CONTINUE 
IS*  RETURN 
EH* 
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Table  1 


RECORD  STRUCTURE  OF  LANDSAT  MSS  CCT 


Data  record 


Number  of  bytes 


JSC  header 
Tape  marker 
Landsat  header 

Geometric  transformation  record 

Calibration  tables  for  band  4 
n  n  n  ti  ^ 


If 

II 

IV 

ft 

II 

II 

If 

tl 

If 

II 

II 

IV 

It 

II 

If 

II 

Tape  marker 
Data  record  1 
Data  record  9144 
Tape  marker 


Table  2 


THE  DATA  TAPES 


RAE  No. 

Satellite  date 

Scene 

Process  date 

Type 

002299 

002272 

001785 

000183 

001855 

9  July  1977 

9  July  1977 

2  July  1977 

2  June  1976 

6  Nov  1975 

208/18 

208/18 

219/24 

218/24 

218/24 

5  Apr  1979 

16  Sept  1979 

9  Sept  1979 

2  Jan  1978 

30  Oct  1979 

Raw 

New  cor. 
Old  cor. 
Old  cor. 
New  cor. 

ITEMS  IN  JSC 


RECORD 


Byte  Nos. 

Type 

Contents 

2738-2745 

2746-2753 

2885-2886 

2887-2894 

EBCDIC 

EBCDIC 

Binary 

EBCDIC 

Sun  elevation,  milliradians 

Sun  azimuth,  milliradians 
Earth  rotation,  milliradians 
Satellite  altitude,  metres 
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Format  specification  for  Landsat  MSS  system 
corrected  computer  compatible  tape  (1976) 

Format  specification  for  Landsat  MSS  system 
corrected  computer  compatible  tapes  produced  at 
Fucino  (Italy)  (1979) 
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